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ABSTRACT 


A thin screen scintillation theory for a spherical wave is 
presented imder the "quasi-optical" approximation. We cal- 
culate the "scattering angle", the "observed angle", the in- 
tensity correlation function and the temporal pulse broadening 
for the random wave. It is found that as the wave propagates 
outward away from the phase screen, the correlation scale of 
the intensity fluctuation increases linearly while the "observed 
angle" decreases linearly. The calculations are carried out 
for both Gaussian and power-law spectra of the turbulent 
medium. 
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THEORY OF THIN SCREEN SCINTILLATIONS 
FOR A SPHERICAL WAVE 


I. INTRODUCTION 

^enomena such as the twinkling of starli^t, the ionspheric, interplanetary 
and interstellar radio wave scintillations involve the scattering of an electro- 
magnetic wave by a layer of random medium. When the thickness of the irr^i- 
larity layer is much smaller than the distance between an observer and the 
irregularity layer, "thin phase-changing screen approximation" (Salpeter 1967, 
Lee and Jokipii 1975 b, c) can be used to analyze the observed scintillati<ni data. 
Applications of the "thin screen scintillation" theory have been successful in 
interplanetary scintillations (e.g. Cohen et al. 1967), in interstellar scintilla- 
tions (e.g. Salpeter 1969, Lee and Jokipii 1976), and in scintillations of radio 
waves during occultation by comet's tail (Lee 1976). 

The existing theory of "thin screen scintillation" assumes an initial plane wave 
impinging on the irregularity layer (Mercier 1962, Salpeter 1967, Lovelace 
1970, Lee and Jokipii 1975b, c). However, most stars or radio sources are 
compact and emit spherical waves. The "plane wave" assumption applies only 
to cases where (rj/rj) ~ 1. Here rj is the distance between the source and 
the irregularity layer, and r 2 is the distance between the source and an observer. 
For interplanetary and comet's tail scintillations, rj/r 2 1 is satisfied and the 
plane-wave thin screen scintillation theory can be applied. For interstellar 
scintillations, r^rj ~ 0.5 and the plane-wave theory can still be applied with- 
out much risk. However, recent observations of angxUar broadening and tem- 
poral broadening of the Crab Nebula pulsar's radiation indicate the existence of 
a scatterii^ region of plasma turbulence within the Crab Nebula (Rankin and 
Counselman 1973, Lyne and Thome 1975, Vandenberg 1976). In this case, r^/rj 
« 1 and the plane wave approximation can not be applied. In o3*der to analyze 
the observed data correctly, a theory of thin screen scintillation for a spherical 
wave is needed. 

It is the purpose of the present paper to present a thin screen scintillation 
theory for a spherical wave. In Section II, we formulate the scintillation problem 
for a spherical wave. In Section in, we derive the scattering angle and the 
obseiwed apparent angle 6^ of the radiowave caused by the plasma irregularity 
layer. In Section IV, we consider the intensity fluctuations for both weak and 
strong scintillations. In Section V, we develop a theory for the temporal broad- 
ening of pulses. In Section VI, we summarize the results and compare them 
with the plane wave case. The theory developed here will be applied to analyze 
the temporal and angular broadening data of the Crab Nebula pulsar in a future 
paper. 
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n. GENERAL CONSIDERATION FOR A SPHERICAL WAVE 


We consider a spherical wave, of wavenumber k = 2tt/x jukI frequency (a;/ 2rr), 
propagating outwaxti tvom a radio source situated at r = 0, 

= F(r,k)e-‘“* = u(r,k) e-‘ >/r (1) 

Let r = (r, ^ be the spherical coordinates in space. As shown in Figure (1), 
the wave propagates freely until it hits the "thin phase- changing screen" situ- 
ated at r = r j. After passing the screen, the phase of the wave is randomized 
and is characterized by the function (t^, 0 ). Then we have 

u(r - d,4>,k) = (2) 

where we have normalized u (r, k) such that 

u(r<rj-D, k) = 1. 

We will consider in this paper the propagation of a wave with o> » ojp , the 
plasma frequency of the plasma medium, thrcmgh the plasma layer. Thmii the 
phase function 4^ is related to Hie electron density fluctuation 6 N^ ( r) 

by 


\( 0 , 4 >) = (-47Tr^/k) p 5N^(r*. 0,(^)dr' (3) 

.'r j-D 

where is the classical electron radius and D is the thickness of the plasma 
layer. ]b deriving Equation (3), we have n^lected the diffraction effect inside 
the plasma layer between r^ - D and r j (Lee and Jokipii 1975b). 

In order to proceed further, we specify the two-point correlation function 



of the fluctuating electron density §N^ ( r), from uliich the statistical proper- 
ties of the phase function ( ', 0 ) can be"^ obtained through Equation (3). As- 
sume that the plasma medium is statistically homogeneous and isotropic. Ilien 
the correlation function ( r) depends o^y on r = I r |. It is convenient here 
to work with the spatial power spectrum ^(q) of the fluctuating electron density, 
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which is related to the correlation function through 

Pn (d) = Pn ( r ) exp (i q • r ) . (4) 

We consider two types of power spectrum for (q): a Gaussian spectrum atid 
a power-law spectrum. For a Gaussian spectrum, we write 

Pn (d) ” Bq (" / 2) (Gaussian) (5a) 

where a is the correlation scale of the fluctuations in electron density. For a 
power-law spectrum, (q) has the following form: 

Pw(q) ' Bp ( 1 + ^ cxp(- / 2 ) (Power- law) (5b) 

where L is the correlation scale of electron density fluctuations and f is the 
cutoff (or inner) scale. Usually f « L, and Equation (5b) for a = 11/3 is the 
Kolmogorov turbulence spectrum. Note that in each case B may be related to 
the mean square electron density fluctuations. We have 

Bq - <5N^2> (ga) 

Bp = <b > r (a '' 2) ' T(a / 2-3 2) (6b) 

where < 6 N/ > is the rms fluctuation electron density. 

Define ^ = (0, ^), ^ = 1 = (y^ + <f>^)^. From Equation (3), we obtain the 

correlation of phase fluctuation, 

= (4.- r, k)^D • (7) 

where 
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f 




1 r* 

2j ( *" 1' ^ 1 ' ® Ng ( r j + r . (K ’Pi ■*■ 0 )^dr ( 8 ) 

— OB 

Here = (in r^/k)^ A,^(0) is the rms phase fluctuation and is the normalized 
phase function, i.e. (0) = 1. In obtaining Equation (7), we have assumed 
that D is much greater than the correlation scale a or L. Assuming D « r p 
we then have from Equations (5a), (6a) and (8), 


Bg 

An(?) “ ; : exp(-r2,f2/2a2) (9) 

in a* 


for the Gaussian spectrum. Now conside.' the power-law spectnun in Equa- 
tion (5b). We have for r , » r, by neglecting the effect of the cutoff at 

q > (- -i, 


) 


Bp(rt -; Lf K^(r,^ L) 
22*^ --L^rOi + 1) 


( 10 ) 


where u = ^ " 1» 2>x ^ 0, and K,^ denotes a modified Bessel function of 

the second kind. One can further show that for L » r, - » f , 


An(^) - 


Bp 


ii, (a - 2) L2 


1 - 


'■(i) 


However, for r, - < f. 


A^(.f) - An( 0 ) 


An'(O) , A^^>(0) 


-.2 4 


4! 


+ 


where 


( 11 ) 


( 12 ) 

I 


4 


r 


I 


j 


4 


4 




An(0) = 


3Bp 
20/r L2 ’ 


a;(0) 






t2 


^•0 


1!p 

16 tr 


L-n^3(f2/2)-‘^® r,2. 


HAd 


A4^>(0) 


d^A^(-) 


d5" 


33p 
64 7T 


L-11/3 (£2/2)-’^6 r0 r|». 


<f=0 


for a = 11/3. 


We shall be interested in the statistical properties of the fluctuating wave at 
r = rj. After passing the screen, the radio wave propagates in free space and 
the wave function F (r, k) satisfies the following equation, 


(V2 + k^) F(r, k) - 0. 


(13) 


Equation (13) can also be written in terms of u as 


2ik 


u 

r 



1 

r ^ sin & 




1 u 


0. (14) 


We are interested in situations where fluctuation scales of u in the radial direc 
tion are much larger than wavelength of the radiation \, so that 
Equation (14) can be neglected, x'he omission of the term u/ dr ^ is called 
the "quasi-optical" or "parabolic" approximation (Lee and Jokipii 1975a). As- 
suming the scattering angle caused by the random phase screen is small. Equa- 
tion (14) can be simplified as 


2ik 


d U 


B r 


+ 




0 


(15) 
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for waves propagating in the direction & -n/ 2, -p -p: 0, Let the observer be 

situated near (r = rj. >■> = t?/2, ^ - 0). We then have from Equation (15), 


u(rj, ■ 




d0, c 




(16) 


where s = rt r 2 /(r 2 - r,). Equation (16) relates the wave function u at r 2 to 
the random phase function <t>^(^\ 7) at r,. We will use Equation (16) to calculate 
the scattering angle, the intensity correlation, and the pulse broadening in the 
following sections. 


m. ANGULAR SCATTERING 

Consider the Fourier transform of the random wave function u(rj,>f ) observed 

atrj, 


u(r2.B) 


(k r2>^ 

4 T7 ^ 



1 2 ) exp ( - i k r 2 0 ‘ 



(17) 


It follows that 


* 


u(r2. 



dB u(r 2 , B) exp(i k f 2 B 


■f)- 


( 18 ) 


For small scattering angles, it can be shown from Equations (11) and (18) that 
u(r 2 , 0) is the Fourier component of the random wave u, propagating in 7!.e 
direction with angle B relative to the radial direction r. We call uirj.B) the 
"angular spectrum" for a spherical random wave. The plane-wave case has 
been treated by Booker, Ratcliffe and Shinn (1950). 

I«t the tvo-point correlation of the random wave u (r j, " 2 ) ‘ 

I2) >2. f 2 5)^ 


where * denotes comple.x conjugate. It may be shown that for the random wave 
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u, the probability of distribution of intensity over angle is related 

to Tj by the relationship 


'Hr 2- ?) 




d-j rj(rj. r) exp(- i • ft) 


( 20 ) 


(Booker, Ratcliffe and Shinn 1950, Lee and Jokipii 1975a). ^{V 2 t^) is the "angu- 
lar power spectrum" of the random wave u. £qu£d.ion (20) states that "the corre- 
lation function of a random wave is proportional to the Fourier transform of its 
angular power spectrum," \n^iich is called the Wiener's theorem (Wiener 1930). 

From Equation (16), we obtain 



\»here s = TjTj /(rj " r,). in deriving Equation (21), we have assumed the 
Gaussian statistics of the random phase function or equivalently, of the 
electron density fluctuation :■ N^, It can be shown from Equation (21) that 



Let the characteristic scale of ! be For £ 1, we have from Equa- 
tion (22), 


■'» r I • for Gaussian spectrum 


- L r j , for power-law spectrum. 


(23a) 
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For > 1, we define such that 


[l“P^(5c)] ^ 1. 


Then we have for a Gaussian spectrum 


VI a 


(23b) 


aivl for a power-law spectrum 



> f. 


(23c) 


gL'^-2(j»2 2)^ 2 


(- 2)1 ( 2 - 


if r, < f- 


(23d) 


The angular power spectrum , (r,,0) at r, is also related to the correlation 
' • :) ^ equation similar to ^uation (20). From Equation (22), we have 


; (f2. 0) 



(24) 


Let the characteristic angles at r^ arid r^ be and 0^^ respectively. Then 
we have 
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II ! i 


e 


0 





(25) 


In Equation (25), we identify as the "observed angle" at and as 
the "scattering angle" due to the plasma Irregularity layer situated at r,. 

Note that the observed angle is smaller than the scattering ruigle by a factor 
of (rj / r^). One would obtain the relationship between B^ and in Equation (2b) 
from the geometry as shown in Figure (1) under the assumption of the validity of 
geometric optics. Note that for a plane-wave, b„ and B^ are of the same value. 

One can also calculate the mean square angles which are defined as 




(26) 


for I = 1, 2. Note that B = We then have 

For a Gaussian spectrum, we have 


(27) 


<«>€> - <^.S> ; K 


I , 


(28a) 


For a power-law spectrum w'ith - 4, L - ’ C , we have 


<«/> 


- ,2 

-) ,0 


2k = 


L' d'2 




(28b) 


IV. THE FLUCTUATIONS IN INTENSH Y 

In order to calculate the intensity correlation function of the random wave u of 
wavenumber k, we consider the fourth moment I ^ which is defined as 
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1 


I 


1 


I 


fZ^ X 'V. 'V 


(29) 


The initial value of T’^ at r^ can be obtained from Equation (2) as 


l4(rt.^a.^>) = {-^0 pP0(O) - 2P^(£J - 2P^(|^,) 




(30) 


+ + ^/?) ^ P;^( "a “ - /O] } 

'X., 'X. 'V. J ^ 


The intensity correlation Pj is related to by 


P,(r.c) - r4(«-.p. 


b. 


0 ) - 1 . 


(31) 


The scintiilation index is defined as 


m 2 -- P,(r,0) 


(32) 


Assuming that the plasma irregularity layer is statistically homogeneous, we 
have from Equation (16) 




— • V, r. 

k f2 ■'-P^ 


(33) 


for r ^ Tj* Define 


• "7, 


" r ■ V 


Equation (33) can also be written as 


■ 4<P' b «) 




, r, 


- V, • Vf P4 . 


(33') 
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\\iiich is Equation (65) of Lee and Joklpii (1975c) for the plane-wave case if t, 
is replaced by z. 

From Equation (33) or (33*), we have 


^"4 (^2» 


k 

2v (ri^ 



1^4 ( t / i . 


exp 


1 k 


(V2 ~ ^1) 




;) 


Let Pj (r, q) and (r, q) be respectively the Fourier transforms of Pj 
and Pj,(^)ln the transverse coordinate c. We have 

Pi(r,q) " (2 /0'^JJPl(r.O d£ 


P^(r.q) ^ (2v,)-'JJp^(£)e-‘^'^ d£. 


It can be shown from Equation (34) that 


Pi (r2- q) 


4 /o^P^(q) sin^ 


(‘‘2 " ri) q^ 
2k fj f2 


for weak scintillations where << 1 (Salpeter 1967, Jokipii 1970). 

For small scattering angle, the local transverse Cartesian coordinate 
lated to .f approximately as; p - r Let the Fourier transforms in 
and P^ be respectively. 


(34) 

Is a . 

(r. ^') 


(35) 


p is re- 
p of P, 
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and 


%(0 



Ap 


Then Equation (35) can be written as 


QiCtj.k) - 4 < p ^ 


^ sin^ 


ri(f 2 - r,) (« fj/r,)* 
2 k r~ 


QiCr, ^/r.) 


(35') 


where x = U | . From Equation (35'), we see that the scale of Fresnel zone for 
a spherical wave is 


pf 



(36) 


(c.f. Born and Wolf 1959, page 372). Consider two limiting cases for p^: 
(i) For Tj / r^ 1, we have 


Pf y. [" (rj - Ti)] 

Here (rj - r j) is the distance between the phase-changing screen and the ob- 
server. The plane-wave approximation is applicable in this case. 

(ii) For rj « rj, we have 


Pk - (37b) 

Here r j is the distance between the radio source and the screen and the plane- 
wave approximation is not applicable in this case. 
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Let pj be the characteristic scale of P, (r,^). From Equation (35’), we have 
(i) For Gaussian spectrum. 


pj rja/rj; (38a) 

(ii) For power-law spectrum. 


Pi - 

if 

L > Pf > f 

(38b) 

= rjf/ri . 

if 

e > Pf- 

(38c) 

.1 

r 

, if 

Pf > L. 

(38d) 


Define rjj and the characteristic distance r^ respectively as 

*■12 “ ri(r2 - ri)/r2 (39) 

and 

= k /0Q (Gaussian spectrum) (40a) 

= (power-law spectrum, fj > f) (40b) 

= kL“-2 (power-law, ri^^<£). (40c) 

It can be shown from Equation (35') that when 

(i) 00 ^ 1 (^^^) 

or 

(ii) ‘/’o ^ ‘■12 fc' (^^^) 
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the scintillation is weak « 1) and Equation (35') gives us the power spec- 
trum of intensity fluctuations. 

For > 1» r ,2 » r^, one obtains from Equation (33) that 


Pi(r2*^) • 

[- 200* [l - P^(£)]} - exp [- 200^] 

(42) 

(Lee and Jokipii 1975c). In this case, the scintillation is strong (M 
scale of intensity fluctuations can be obtained from Equation (42) 

1). The 

Pj = 2"^^k"*0Q'* = i^/%) • (tj/tj) (Gaussian) 

(43a) 

^ 0-1 

(power-law, L > r j > 

f) (43b) 

^ 2-*'^ k“* 0-‘ 

(power-law, < f) 

(43c) 


Mtiiere 0^ and are given by Equations (23) and (25). From Equations (38) and 
(43) , we find that pj « (rj/rj). Thus the intensity correlation scale expands 
as the random wave propagates outward spherically. For the region where 
rj 2 r^, the correlation scale Pj changes smoothly from the values given by 
Equations (38a, b,c) to those given by Equations (43a, b,c) (Lee 1974). 


V. TEMPORAL BROADENING OF PULSES 

The temporal broadening of pulses has been considered by Lee and Jokipii 
(1975b) for the plane-wave case. We consider here the pulse broadening for a 
spherical wave. Let R(r, t) be the probability distribution of pulse intensity as 
a function of time t measured at r. Define the second moment of the random 
wave u between different wave-numbers as 

r(r.|, Ak) </u(r, ?J. k j) u*(r, k 2 ^ (44) 


where f f , and Ak = kj - k^. It can be shown that 
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R(r, t) 


<45) 



J dAk r(r, 0, Ak) exp (- i 6k ct) 


(Lee and Jokipii 1975b). A transport equation for F (r» Ak) can be obtained 
from Equation (IS), 


3F(r,^, Ak) 
+ 


V/r*0(k») 

2k2r* I 


0 . 


(46) 


We assiune Ak is small and | Akj « k and neglect terms of order Ak^ in 
Equation (46). The initial value of F at r ^ can be obtained from Equation (2) as 


F(rp£,Ak) - exp|- ^ <^0^1 exp ^.^2 - P^(£)]} . 


(47) 


By employing the same technique used by Lee and Jokipii (1975b), we obtain 
from Equations (45, 46, 47), 


R(f2. t) ^ R,(t)WR2(r2, t) (48) 

where denotes convolution of two functions in the variable t, and 

Rj ( t ) " (c k /Po ) tt" fxp (- k 2 1 2 /0^) (49a) 


R2(t) 


2t 7 c r j 
^ 2^^ 2 ” l ) 


r^.e 


' 2cr,t - |\ 
/ 2(«’2 - UU )' 


t > 0, 


(49b) 


= 0 


, t < 0. (49b) 


In Equation (49b), P (rj.B) is the probability distribution of wave intensity over 
angle as discussed in Section III. 
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K 1 (t) is the broadenii^ due to the "pure refraction" effect because it gives the 
effect of the differing transit times of different rays due to the varying iiaiex of 
refraction. Note that R ^ (t) is a symmetric Gaussian function with characteristic 
time 


tj = 0(,/ck. (50) 

Rj (rj , t) gives the pulse profile due to the "pure diffraction" effect. From 
Equation (49b), we obtain the characteristic time scale tj of R 2 as 


f2(r2 ~ ri) 

2 c r. 


02 


~ «-i) 
2 c r. 




(51) 


where 6^ and 0^ are respectively the "observed angle" and "scattering angle" 
as given by Equation (25). The forms of Rj (t) for a Gaussian irregularity 
spectrum and a Kolmogorov power-law spectrum are exactly the same as for 
the plane-wave case and are given by Figure (6) of the Lee and Jokipii (1975b). 
For the Gaussian spectmm, R 2 (f) is an exponential function. But for a power- 
law spectrum, R 2 (t) does not have an exponential shape. The fact that R 2 (t) = 0 
for t < 0 indicates that all the scattered rays are delayed. 


VI. DISCUSSIONS 

From the results obtained in the previous sections, we find that the plane- wave 
thin screen approximation is no longer applicable for cases where r j / Tj « 1. 
We summarize the thin screen scintillation results for a spherical wave as 
follows. 

(a) There are two characteristic angles, namely, the scattering angle 0g 

and the observed angle 0^ . The scattering angle 0g is determined by 
the amount of irregularity of the plasma layer and does not depeixl on 
rj or r^. However, the observed angle 0^ decreases as the wave 
propagates away from the scattering screen and is equal to (r, / r^). 

In the plane wave case, 0,, = Bg* The analysis of scattering data of the 
Crab Nebula Pulsar by Vandenberg (1976) is incorrect since the author 
uses 0Q as the angle related to the amount of irr^ularity in the plasma 
layer. 

14 

(b) The scale of Fresnel size ,< > is [ rj (r 2 - Tj ) >^ / r 2 1 instead of 
l(r 2 - r,) ‘ as in the plane-wave case. 


16 


I 


(c) The Intensity correlation scale p, is proporticnal to (r^ / r^ ) and 
increases as the wave propagates outward away from the screen. La 
the plane-wave case, pj reaches a steady value for strong scintilla- 
tions. 

(d) The conditions for strong scintillations are > 1 and r^j « ( 1^2 ' 

rjl/rj > r^. For plane-wave case, the conditions are: > 1, z * 

(rj - r,) > r^. 

(e) The characteristic time scale tj of pulse broadening due to the diffrac- 
tion effect is given by 


2 c r , 


0n 


~ «~l) 
2 c r, 




Finally we note that an observer usually measures the correlation time r, of 
intensity fluctuations. Let the characteristic scale of P, in Equation (42) be ^ , . 
Then we have j = pj / Tj. Asstune the fluctuations in the plasma layer are 
frozen in and the plasma layer moves at an angular velocity i = d^ / dt trans- 
verse to the line of sight to the pulsar. Then 


T 


I 



(•"2 ^ • 


(52) 


If I is caused by the motion of the plasma layer with a transverse velocity v j 
relative to the line of sight to the pulsar, then cj = v, /tj , And we have 


T 


I 


‘‘l 


V, rj 


(53) 


If f is caused by the transverse motion of an observer with velocity v^ , then 
i = Vj/rj. We have 


T 


I 




' 2 - 


(54) 


For example, if the temporal broadening of pulses from the Crab Nebula pulsar 
is due to Filament 159 as proposed by Vandenberg (1976), then rj / rj = 10"^ , 
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V) == 30 km/sec and p, s 2 x cm. One obtaiiw from E<|uation (S3), = 

0.67 second, comparing to Tj = p, /v, - 2 hours as would be obtained from 
the Incorrect plane-wave theory. 

In a later p^er, we will apply the tiieory developed here to analyze the tem- 
poral and angular broadening data for the Crab Nebula pulsar. 
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FIGURE CAPTION 


Figure 1. A schematic sketch of the thin screen diffraction problem for a 

spherical wave. The random medium is confined to a thin layer of 
thickness D (from r = r^ - D to r » rj). The spherical wave 
e* (*''■"*)/ r propagating from the source (r = 0) hits the "thin 
screen". After passing tlie screen, the phase of the wave is ran> 
domized and is characterized by the function (-> , <P). A ray path 
SAO is shown to demonstrate that 0^ = r^ 6g / r^. 
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